


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1982 


Thermal calibration of satellite infrared 
images and correlation with sea-surface 
nutrient distribution. 


Pereira e Silva, Vitor Martinho F. 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/20377 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sha Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


| 2. 


i) 
? 
R ' i 
' 
any 
if 1, 
sf 
i 
sf ' 





DUDLEY KMOX LIBRARY 
AVAL POS” ARmUETE SCHOOL 
ONTEREY, CALIF, 93949 

















NAVAL POSTGRADUATE SCHOOL 


Monterey, Galifornia 





THESIS 


THERMAL CALIBRATION OF SATELLITE 
INFRARED IMAGES AND CORRELATION WITH 
SEA-SURFACE NUTRIENT DISTRIBUTION 


by 


Vitor Martinho F. Pereira e Silva 


June 1982 





Thesis Advisor: E. Traganza 


Approved for public release; distribution unlimited. 


7205760 


| 
| 








UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


REPORT DOCUMENTATION PAGE 


2. GOVT ACCESSION NO. 





READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


3. RECIPIENT'S CATALOG NUMBER 


5. TYPE OF REPORT &@ PERIOO COVEREO 
Master's thesis; 
June..1982 


6. PERFORMING ORG. REPORT NUMBER 


| @ CONTRA OR GRant NUMBER( 2) 










4. TITLE (and Subtitie) 


Thermal Calibration of Satellite 
Infrared Images and Correlation with 
Sea-Surface Nutrient Distribution 













AUTHOR 8) 












Maaeer MartimhoweF. Perera e Silva 





0. PROGRAM ELEMENT. PMOJECT TA 
AREA & WORK UNIT NUM@ERS | = 


12. REPORT CATE 
June 1982 
13. NUMBER OF PAGES | 
64 

. MONITORING AGENCY NAME & ADORESS(1I1 different irom Conitroiling Office) 18. SECURITY CLASS. (ol (hie report) 
Unclassified 





9. PERFORMING ORGANIZATION NAME ANO ADORESS 


Naval Postgraduate School 
Monterey, California 93940 








CONTROLLING OF FICE NAME ANDO AOCORESS 





V1, 


Naval Postgraduate School 
Menterey, California 93940 









Se. OECLASSIFICATION/ DOWNGRAOING 
SCHEOULE 


16. OISTRIBUTION STATEMENT (of tate Repert) 


Approved for public release; distribution unlimited. 


17. OISTRIBUTION STATEMENT (of the edetract entered in Block 20, fi ditterent trem Report) 


18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue om reverece side ii neceesary and identify wy bieck rumober) 


satellite infrared imagery; radiative transfer theory; 
sea surface temperature; nutrients; nitrate; phosphate; 
upwelling; mapping 


20. ABSTRACT (Continue an reverse side if nececeary and identify by bieck maeber) 


Satellite infrared imagery off the California coast, near 
Pt. Sur, show thermal patterns associated with an upwelling 
center; the patterns frequently curl cyclonically when inter- 
acting with the warmer California Current. This pattern shows 
sharp thermal fronts, easily identified in satellite IR images, 
that are strongly correlated with nutrient fronts during the 
early stages of upwelling. With sea truth data available, it 





DO en's; 1473 cortion oF | Nov 6818 OBSOLETE 
S/N 0102-014+ 6601 | 


NCL PILED 


SECURITY CLASSIFICATION OF THIS PAGE (Prem Deore Entered) 





UNCLASSIFIED 


SOCUINTY CL ABSIFIC ATION OF Tmih PAGE Wrem Nore Entered. 


#20 ~- ABSTRACT - (CONTINUED) 


was feasible to calibrate satellite derived sea surface 
temperature, by applying radiative transfer theory, and to 
infer nutrient concentrations from their linear inverse 
correlations with temperature. Thus, 1t was possible to 
calibrate satellite thermal fields to produce maps of 
nutrient distributions. When the inferred relationships 
were applied over representative regions of the upwelling 
center, standard deviations of 0.5°C, 1.7 uM and 0.1 uM 
were computed for temperature, nitrate and phosphate, 


respectively. 





SS lle UNCLASSIFIED 
oN 01022914-6601 SE CYAITY CLABGBIPIC ATION OF Twis Pad (Oren Sara Insered) 


to 











Approved for public release; distribution unlimited. 


Thermal Calibration of Satellite 
Infrared Images and Correlation with 
Sea-Surface Nutrient Distribution 


by 
Viseaeicdtaeimine F. Pereira e Silva 


Lieutenant, Portuguese Navy 
B.A., Portuguese Naval Academy, 1976 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN OCEANOGRAPHY 


from the 


NAVAL POSTGRADUATE SCHOOL 
June 198z 





ABSTRACT 


Satellite infrared imagery off the California coast, 
near Pt. Sur, show thermal patterns associated with an up- 
welling center; the patterns frequently curl cyclonically 
when interacting with the warmer California Current. This 
pattern shows sharp thermal fronts, easily identified in 
Satellite IR images, that are strongly correlated with nutri- 
ent fronts during the early stages of upwelling. With sea 
truth data available, it was feasible to calibrate satellite 
derived sea surface temperature, by applying radiative trans- 
fer theory, and to infer nutrient concentrations from their 
linear inverse correlations with temperature. Thus, it was 
possible to calibrate satellite thermal fields to produce 
maps of nutrient distributions. When the inferred relation- 
ships were applied over representative regions of the upwelling 
center, standard deviations of 0.5°C, 1.7 uM and 0.1 uM were 


computed for temperature, nitrate and phosphate, respectively. 
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I. I[NGRODUCTION 


This thesis is a study of the possibility of calibrating 
ied =red Satelite images from limited in-situ temperature 
data, and the possibility af whether chemical (specifically 
nitrate and phosphate concentrations) mesoscale variability 
may be inferred from corresponding satellite detected thermal 
patterns in an upwelling zone [Traganza, 1980]. 

Very high resolution infrared radiometry, with a spatial 
and amplitude resolution ~1.0 km and ~-0.5 °C, respectively, 
have been successfully used in locating and sensing sea 
surface temperature contrasts in upwelling zones [Traganza 
pemol., L980; Bowman et al., 1977]. Due to the development 
in radiative transfer theory and satellite imagery it is 
possible to calibrate: this satellite derived temperature 
when sea truth data are available, and the atmospheric mois- 
fume Cdis@rmeution is uniform. 

Surface temperature and nutrient maps can be generated 
using these corrected satellite derived temperatures, and 
nutrient to temperature correlations, because active upwelling 
systems in an early stage of development are expected to have 
strong inverse linear correlations between nutrients and 
temperature. The recurrent formation of "a cyclonic upwelling 
system" off Pt. Sur, California [Traganza, Conrad and Breaker, 
1980] offers a unique opportunity to simultaneously investi- 


gate the accuracy of satellite generated surface maps for 


PL 





temperature and nutrients, and the influence that coastal 
upwelling has on surface nutrient concentrations when it 
brings nutrient-rich waters into the euphotic zone. This 
appears to be a major factor regulating the standing stock 
Bm@sprocuction of phytoplankton in coastal waters [Eppley et 
Bue, 42/9). And, because phytoplankton are the primary pro- 
ducers in the pelagic ecosystem, supporting the foodweb, it 
will be of interest to obtain a wider and instantaneous pic- 
ture of the nutrient concentrations, which only a satellite 
Bam give, although indirectly. 

Thermal gradients detected by IR imagery identify potential 
Site of high biomass which may be confirmed by Coastal Zone 
Color Scanner (CZCS) data [Traganza, 1981]. This data can 
reveal subtle variations in the concentration of phytoplank- 
ton pigments and has a potential application for the study 
of large-scale patchiness in phytoplankton distributions 


Mievis et al., 1980]. 





ii MEREODS 


mee CRUISES 

Sea truth data were obtained from the cruises conducted 
aboard the Naval Postgraduate School's research vessel R/V 
ACANIA on 10 and ll June 1980 and 29 October 1980. 

To ensure coverage of the cyclonic moving feature updated 
satellite information regarding its approximate location, size, 
center and orientation, was sent prior to each cruise and 
sometimes during the cruise by Mr. Breaker, the staff ocean- 
Ographer at the National Environmental Satellite Service (NESS) 
Station, in Redwood City, California. R/V ACANIA positions 
within this "cyclonic upwelling system," centered a few 
kilometers southwest of Point Sur, were determined by LORAN-C 
which has good coverage in this region. Temperature and nutri- 
Pieomvere Measured in Situ every two minutes at a depth of 
2.5 m. Seawater temperature was sensed by a thermistor, and 
recorded continuously on a strip chart. For calibration and 
monitoring of the equipment sea surface temperature was also 
measured by bucket thermometer. Surface nutrient concentra- 
tions were recorded and analyzed every two minutes according 


to the Technicon Industrial Method [Hanson, 1980]. 


B. SATELLITE IMAGERY 
Three satellite images of the region offshore Central 


California were analyzed and calibrated. They were chosen 


ee 





because they were cloud-free over the upwelling feature and 
taken during or very close to the cruise time. 

The infrared (10.5 to 12.5 uM) data recorded in the 
satellite orbits were analyzed at the NASA-Ames Research 
Center, in a computer equipped with the Interactive Digital 
Image Manipulation System (IDIMS), which is a comprehensive 
software package developed by Electromagnetic Systems Labora- 
tories, Inc. With this software package, the author could 
"zoom" magnify and pseudo-color the displayed image under 
myOvSstick Control. Some 8=by-10 inch color Polaroid prints 
were taken, as shown in Plates 1, 2 and 3. A computer print- 
out (picprint), with the recorded count values (radiometric 
units of measurement) was also obtained (Figs. 1, 8, 15). To 
navigate on it, and because there were no geographical coor- 
dinates on the picprint, the author used a computer program 
developed by Lundell (1980). The purpose of this program 
is to determine the line and pixel number of a geographical 
location given the locations latitude and longitude, of a 
landmark location (line, pixel, latitude and longitude) and 
the period and inclination of the satellite orbit. 

With the image navigational problem solved, it was possi- 
ble to draw portions of the R/V ACANIA cruise tracks on the 
picprint in order to compare the digitized satellite tempera- 


tures with the lie St values. 
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C. INFRARED RADIOMETRY 
 Teriperatume Degradation 

The usefulness of satellite observations depends on 
the ability of the radiometer to view the sea with little 
error introduced by the atmosphere. This error, which can 
be substantial even ina relatively clear air, 1s mainly in- 
duced by atmospheric attenuation of the infrared radiance 
due to absorption by water vapor. According to Maul and 
Sidran (1973), this effect varies from small values for arc- 
meen aucm@spiercs to Eairly high values, ca. 10°C for tropical 
atmospheres. In mid-latitudes, this attenuation may introduce 
Steers on the order of 2 to 6°C. 

Water, both as vapor and as clouds, is the primary 
source of error. Even in cloud-free atmospheres, where scatter- 
ing of infrared radiance is negligible, the magnitude of ab- 
sorption and the variability of the moisture field presents 
a severe problem. Theoretically, the influence of water vapor 
can be calculated if the absorption coefficient of water vapor 
Me well aeethe vertieal temperature and humidity profiles are 
known. However, this requires radiosonde data for a study 
area, which is rather the exception. Even using an isolated 
sounding the temperature accuracy is limited to 1°C [Maul and 
Sideam, 19)7 3a. 

Failure to detect the increase in humidity, due to 
Clieuge, leag@e to 2 te 5°C errors in the calculations. Even 
in cleareatcmecprenes, the water vapor can reduce the apparent 


sea surface temperature by 4 to 8°C. Furthermore, the 


iS 











presence of moisture markedly attenuates surface temperature 
feadvets deross am oceanic front [Maul et al., 1978]. 

Because the influence of clouds is the most difficult 
variable to estimate and to correct for, the cloud-free satel- 
lite observations, present during this study, minimized the 
atmospheric errors and simplified the radiative transfer 
equation because scattering of infrared radiance could be 
neglected. 

Adding to these atmospheric errors and increasing the 
total deviation between satellite derived temperatues and 
B@esitu Cemperatures, are errors. Originating at the sea sur- 
face, due to the very strong absorption (a) of infrared 
photons on it. This produces a very thin optical depth 


CZ mr ce, wom ene Crdereorsca.,) 0.006mm for the ca. 10.5 


Iya 
uM to 12.5 uM IR spectral window. AS a consequence, the IR 
radiation measured by the radiometer originates in the upper 
0.006 mm of the sea surface, and may not represent slightly 
deeper temperatures, as those obtained by the thermistor at 

2.5 mon the R/V ACANIA cruises. The very thin "skin" at 

the top of the ocean iS in a gquaSi-~laminar viscous sub-layer, 
dominated by molecular fluxes, which are the reason for the 
discrepancy between the sea surface temperatures sensed remotely 
Byettcideeenanc DY ai-situ methods. The difference is usually 
Clete OrdemNOt ice “8eo5to 1°C, with the in-situ values 

warmer [Stewart, 1979]. 


All these errors originating in the atmosphere and at 


the ocean surface will be corrected indirectly in this study 
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by correlation between temperature derived radiances from 
aoeewrite and in-situ data. The author will assume negli- 
gible both water vapor gradient and temperature discrepancy 
(between satellite derived and in-situ temperatures) gradi- 
ent within the area. 
me Planck Pumetion 

The sea surface emits thermal infrared spectral radi- 
ance (L) as a function of the wavelength (A), the sea surface 
temperature (T) and zenith incidence angle (8) according to 


the equation 
pawn, 1) = Seiten) <6 (i, T) 
where 8(A,T) is the Planck Function, 


2 
Ba AGL) Ons SL ieee seaward denen 3 


1>. [exp (he /kTA) -1] 


where 
h = Planck's constant = 6.626196 «10 °* joule.sec 
A = wavelength of emitted radiation, meters 
c = speed of light = 2.997925 x10° m-sec > 
ee ; ” =23 . worl 
k = Boltzmann's constant = 1.380622 x 10 joule*K 
T = temperature, K 


The Planck Function gives the spectral radiance emitted 


by an ideal blackbody. However, the sea surface 1s not an 


t- 
~J] 





ideal blackbody and emits less radiation. The factor ¢«(A,9) 


accounts for this fact, and is called Emissivity, ¢(A,9) 


8(sea surface, K) 


e(A,8) 3(blackbody, K) 


Because its averaged value over the IR spectral window (10.5 
um to 12.5 um) remains near eared £Or ge Came Oe to 40", 
this value was used in the computations of radiances. 
ce RaGiaeivelivanster Theory 
Absolute measurements of sea-surface temperature could, 
in theory, be produced from solutions of the radiative 


transfer equation, 


eee CO mea (DL 


where 


L(T,) radiance measured by the satellite radiometer 
in watts-steradian *~-m™ 

ii = temperature obtained by inverting the Planck's 
function for the measured value of radiance 

in K 


total radiance emitted from the surface 
at temperature Te in watts:-steradian™+-m 


iy (CE =3 


atmospheric transmittance between the 
peeesune levels ca. 0 (top of the atmos- 
phere) and Po (sea level) 


IE = path radiance due to the isotropic thermal 
emission of photons by the atmosphere along 
the propagation path from the sea surface 
to the radiometer in watts-steradian~!-m7 
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using detailed atmospheric observations and multi-channel 
satellite radiometers. In this study, such data were not 
available, and therefore, a simpler scheme was used of 
necessity. To derive the parameters t, the atmospheric trans- 
[ttanee, and L_, the path radiance, both sets of satellite 
Piemaiiesttu temperature observations in the observed frontal 
zones, where thermal differences are large enough to reduce 
the relative errors due to satellite measured thermal noise 
(e.g., caused by variable moisture over the area) and due to 

a coarse AVHRR (Advanced Very High Resolution Radiometer) 
temperature resolution of 0.5 Ton More precisely, linear 
regression analysis was done on both sets of temperature 
derived radiances versus elapsed distance along the cruise 
M@eeex, In frontal regions (Figs. 2, 9, and 16). Assuming that 
the sea surface temperature gradient is much larger than the 
moisture gradient, then |7L,| >> |VL,|, then by the radiative 


transfer equation 


tT = VL/VL5; Te iSecoeek! 
within a very few percent. The ratio of the slopes of the 
regression lines will give then, the transmittance, +t. The 


other atmospheric correction factor in the radiative transfer 


* 

Although the inherent quantization of the AVHRR is ~-0.2 K 
at 285 K, the NESS station in Redwood City can only record 
eight bits of the 10 bits transmitted by the satellite, thus 
degrading the resolution to 0.5 K in the data available for 
EES meucy . 


eo 





equation, the path radiance, L,, can be estimated using the 
Pee CD commen the Yoaxis Of the satellite and in-situ radi- 


ances regression lines, respectively L' and Lo, knowing that 


= ' 
L mx + ZL 


and 


where my ow reo chews hopes OL the satellite and in-situ 


2 
radiance regression lines. Taking in account the relationship 
between the slopes of both lines, t = m,/M5, and substituting 


in the radiative transfer equation, the equation for the path 


radiance is obtained: 


Finally, knowing both transmittance, +t, and path 
radiance, L,, any apparent (measured by the AVERR) sea surface 
temperature was corrected for using the radiative transfer 


equation, with reasonable accuracy (residual std. dev. 0.5°C). 


D. LINEAR REGRESSION AND CORRELATION COMPUTATIONS 

Least squares linear regressions, designed to minimize 
the sum of the squares of the deviations of the data points 
from the straight line of best fit, were performed in nitrate 
versus phosphate, nitrate versus temperature and phosphate 


versus temperature. Regression lines were also generated for 


20 





Eaavences, Computed from in Situ and satellite temperatures, 
versus elapsed distance along the portions of the cruise track 
with strong temperature gradients and nearest in time to the 
Satellite overpass. The slopes and y-intercepts from the 
regression lines were computed from the equation 

UX. Ys = NEX-Y3 


—- aa 


a 
y - i 


be = ey = mx 


The population correlation coefficient was also computed 


for the nutrient regression lines, using the equation 





where oe and ay are the standard deviations of the x and y 


values, respectively. 
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Ttt. ReSvLts 


A. 9 JUNE 1980 SATELLITE IMAGERY 

Several satellite images off Central California prior to 
and including 9 June, showed a "tongue" of upwelling water 
within an early stage of development, near Pt. Sur. 

The best cloud-free image was taken by the satellite NOAA- 
6, in its orbit number 4954, at 0250 GMT 9 June 1980. The 
Meeeaeea (L025 pm to 12.5 um) data recorded in this orbit was 
analyzed, at the NASA-Ames Research Center, in a computer 
equipped with the IDIMS. Plate 1 shows the black and white 
copy of the pseudo-color image processed. A computer print- 
out (picprint) with the recorded count values was also obtained. 
With the identification of several geographical locations 
(Pt. Sur, Pfeiffer Pt. and Lopez Pt.) and using the computer 
program developed by Lundell (1980), the portion of the R/V 
ACANIA cruise track nearest in time to the satellite overpass 
was drawn on it (Fig. 1). Temperatures were digitized along 
eo] Cremer de« and COmMpared with the In-situ values (Fig. 2). 
Beg GOS cumernes Were COmpured for both im-situ and satellite 
temperature derived radiances, in the regions with strong 
Gradienes, Viz., between elapsed distance ca. 4 km to 16 kn, 
30 km to 60 km and 64 km to 83 km. From the comparison between 
Piese limes, avemagec values for the tramsmittance ca. 0.693 
and the path radiance Ca. 2,150,000 Wetec e oS acsol] aie were 


computed. Using these values in the radiative transfer 
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Plawe 1. 





P= oe 


PT, LOPEZ 
BS 3 


NOAA-6 satellite image of the California coast 
for 9 June 1980. Letters associated with color 
grades represent averaged temperatures, over a 
pixel, measured by the radiometer (see Table I 
for radiometric count, temperature, nitrate and 
phosphate values). 
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ipeeeUeremperatusze (Straight line) and 
satellite temperature (line with circles) 
versus elapsed distance along a portion of 
the cruise track of the 9-10 June 1980 
cruise. 
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equation it was possible to correct the satellite temperatures 
associated with the count values (Table I). 

The sea surface nutrient distribution was obtained from 
Pmeo@ercl ation with the in-situ thermal distribution. Linear 
regression analysis between nitrate and temperature, phos- 
phate and temperature, and nitrate to phosphate yielded slopes 
@e —-35.55, -0.27 and 12.79 respectively, and y-axis intercepts 
of 49.27 uM, 4.72 uM and -12.25 uM respectively (Figs. 3, 4 
and 5). The correlation coefficients were r = -0.96, r = -0.96 
and r = 0.98 respectively. From this analysis a nutrient 
concentration was computed for each calibrated satellite 
temperature (Table lL). 

Sea surface temperature nitrate and phosphate maps were 
generated by “zoom transfering" the thermal pattern in the 
IR satellite image to a navigational chart. Each isopleth 
is the average of the temperature or nutrient concentration 


on each side (Figs. 6, 7 and 8). 


B. ii JUNE 1980 AND SATELLITE IMAGERY 

A strong "cyclonic upwelling feature" [Traganza, E.D., 
eieat., (1980); Traganza, E.D., Austin, D., (1961)} is shown 
in the image taken by the satellite NOAA-6, in its orbit num- 
ber 4983, at 0350 GMT 11 June 1980. Plate 2 shows the black 


and white copy of the pseudo-color image processed with IDIMS. 


* e 
This was done using zoom transfer scope, an optical 


device used to superimpose one magnified and linearly 
stretched image on another, so that spatial features may be 
M@e@emstersed by tzactng with a pen. 
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Table I 


Calibration Table of the satellite image for 9 June 1980 


(the '*' represents values below limit of detection) 


Color Satellite Corrected Nitrate Phosphate 


Grade Count Temperature (°C) Temperature (°C) NO. (uM) Ore (uM) 


K 103 Le, OC 14.92 x 0.69 
wo. - 8«=—«édz LA $4563 1 x 0.85 
I 105 12.03 13.68 0.73 iL, 10) 
stay 1.83 LLL 

H 106 11.59 13.06 2.93 1.19 
ee 5 4.05 1.28 

G 107 11.15 iA 5.16 1.36 
Lala 6.30 1, as 

F 108 10.70 11.79 7.44 ey: 
11.48 8.56 1.63 

E 109 10.26 flea 6 9.67 1 
10.84 10.82 1.80 

D LEO 9.81 10.51 11.96 1.88 
10.19 eel 1.97 
C ie 9.36 9.87 14.25 2.05 
9.55 PS 2410 2.14 

B ial 8.91 9.22 ion 55 es 
8.90 1770 232 

A Was 8.46 2) 5 Bz 18485 2.40 
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Figure 3. Nitrate versus temperature with the regression 
line for 10 June 1980. 
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Figure 4. Phosphate versus temperature with the regression 
line for 10 June 1980. 
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Nitrate versus phosphate with the regression 
ime for lO Tune 1980. 
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Figure 6. 





PT, LOBOS 


sea surface temperature maps in °C for 9 June 
1980 inferred from satellite IR imagery. 
Contour interval is one radiometric unit of 
measurement (1 count = 0.5°C). 
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Figure 7. 





PT, LOBOS 


Sea surface nitrate map for 9 June 1980, 
generated by correlation with sea surface 
temperature distribution given by IR imagery. 
Contour interval is one radiometric unit of 
measurement (1 count = 0.5°C). Nitrate con- 
centrations are given in uM. 





Figune 8. 
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Sea surface phosphate map for 9 June 1980, 
generated by correlation with sea surface 
temperature distribution given by IR imagery. 
G@re tour interval is one radiometric unit of 
Medow@pcment. (1 count = 0.5°C). Phosphate 
concentrations are given in uM. 
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Plate 2. 





OPEL PBSPE SLPS IOD EBT L ND” BEEN 


NOAA-6 satellite image of the California coast 
for ll June 1980. Letters associated with color 
grades represent averaged temperatures, over a 
pixel, measured by the radiometer (see Table II 
feewracrenetric Count, temperaturc, nitrate and 
phosphate values). 
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The computer printout (picprint) was also navigated, using 

the computer program developed by Lundell (1980), in order 

to draw on it the portion of the R/V ACANIA cruise track 
nearest in time to the satellite overpass (Fig. 9). Tempera- 
tures were digitized along the cruise track and compared with 
mgenin=situ values (Fig.10). Regression lines were computed 
PeeecOraeiaas tty eamd satellite temperature derived radiances, 
Mimemwe) Feq~ons with strong gradients, viz., between elapsed 
peeeermce Ca. 219 km to 223 km, 223 km to 230 km and 240 km 

to 262 km from the comparison between these lines, averaged 
Pees fOr the Cransmittance ca. 0.525 and the path radiance 
Gam 3,320,000 watts-steradian +-m7> were computed. Using 
these values in the radiative transfer equation, a calibration 
table for the satellite temperatures associated with the count 
values was constructed (Table II). 

The sea surface nutrient distribution was obtained from 
ieomeer relation with thewin-situ thermal distribution. Linear 
regression analysis between nitrate and temperature, phosphate 
and temperature, and nitrate to phosphate yielded slopes of 
-4.26, -0.29 and 14.46 respectively, and y-axis intercepts of 
57.67 uM, 4.85 uM and -13.54 uM respectively (Figs. li, 12 
Bewels). eemecorrelation coefficients were r = -0.97, r = -0.96 
and r = 0.99 respectively. From this analysis a nutrient 
concentration was computed for each calibrated satellite 


temperature (Table II). 
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Figure 10. In-situ temperature (straight line) and 


Satellite temperature (line with circles) 
versus elapsed distance along a portion of 
the cruise track of the li June 1980. 
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Calibration Table of the satellite image for ll June 1980 


(The '*' represents values below limit of detection) 


Color Satellite Corrected Nitrate Phqsphate 
Grade Count Temperature (°C) Temperature (°C) NO. (uM) PO 4 (uM) 


18 105 1A. Oe mS 7 * 0.48 
K 106 11.59 14.48 0.68 
3 107 11.15 13.80 * 0.87 
z 108 10.70 13.10 1.82 1.07 
. 17 6 S27 leaelig? 

H 109 10.26 a. 4.76 ee 
12.06 6.26 1.37 

6 110 9.81 ie Te 7.79 1.48 
11.35 9.28 1.58 

F 111 9.36 11.00 10.78 1.68 
| 10.65 ee 7 1.78 

E Lah? oon 10.29 13.80 1.88 
9.94 15.29 1.98 

D 113 8.46 9.58 16.83 2.09 
9.22 18.36 2.19 

C 114 8.00 8.85 19.94 2.30 
8.49 21.48 pe 

B 115 7.54 s7 12 23.05 2.51 
7.76 24.59 Datel 

A 116 7.08 7.39 26.17 2a? 
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Figure ll. Nitrate versus temperature with the regression 


line for 1l June 1980. 
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Figure 12. Phosphate versus temperature with the 
regression line for 11 June 1980. 
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Figure 13. Nitrate versus phosphate with the regression 
Mine tor 1) June 19310. 
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Sea surface temperature, nitrate and phosphate maps were 
generated by "zoom transferring". the thermal pattern in the 
IR satellite image to a navigational chart. Each isopleth 
is the average of the temperature or nutrient concentration 


on each side (Figs. 14, 15 and 16). 


eee 22 OCTOBER 1980 SATELLITE IMAGERY 

The image taken by the satellite NOAA-6, in its orbit 
number 6967, at 1640 GMT 29 October 1980, shows a different 
cold feature, compared with the one from the June 1980 cruise. 
In this case, there is an elongated "tongue" of cold water. 

Plate 3 shows the black and white copy of the pseudo- 
color image processed with IDIMS, in the NASA-Ames Research 
Center. The computer printout was again navigated using the 
computer program developed by Lundell (1980), in order to 
draw on it the portion of the R/V ACANIA cruise track nearest 
in time to the satellite overpass (Fig. 17). Temperatures 
were digitized along the cruise track and compared with the 
mmoerdevaiues (Fig. 18). Regression lines were computed for 
ene in-situ and satellite temperature derived radiances, in 
Paoeeed Omer with Strong gradients, viz., between elapsed diis- 
momee ca. Zi km to 51 km, 31 km to 70 km and 170 km to 191 km. 
From the comparison between these lines, averaged values for 
a= Crememectamce ca. 0.685 and the path radiance ca. 2,059,000 
watts’ steradian /-m7? were computed. Using these values in 


the radiative transfer equation, a calibration table for 


* 
See previous footnote. 
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Figure 14. Sea surface temperature map rie Or bd 


June 1980 inferred from satellite IR imagery. 
Contour interval is one radiometric unit 
of measurement (1 count = 0.5°C). 








Figume 15. 





PT. LOBOS 





Sea surface nitrate map for ll June 1980, 
generated by correlation with sea surface 
temperature distribution given by IR imagery. 
Contour interval is one radiometric unit of 
measurement (1 count = 0.5°C). Nitrate 
concentratins given are in uM. 
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Frgpee 16. 





PT, LOBOS 


Sea surface phosphate map for ll June 1980, 
generated by correlation with sea surface 
temperature distribution given by IR imagery. 
Contour interval is one radiometric unit of 


measurement (1 count = 0.5°C). Phosphate 
concentrations are given in uM. 
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Plate 3. NOAA~6 satellite image of the California coast 
for 29 October 1980. Letters associated with 
color grades represent averaged temperatures, 
Over a pixel, measured by the radiometer 
(see Table III for radiometric count, tempera- 
ture, nitrate and phosphate values). 
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Figure 18. In-situ temperature (straight line) and 





Satelite temperature (line with circles) 
versus elapsed distance along a portion 

of the cruise track of the 29 October 1980 
cruise. 
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the satellite temperatures associated with the count values 
was constructed (Table III). 

The sea surface nutrient distribution was obtained from 
its correlation with the in-situ thermal distribution. Linear 
regression analysis between nitrate and temperature, phosphate 
and temperature, and nitrate and phosphate yielded slopes of 
-1.27, ~.09 and 13.60 respectively, and y-axis intercepts 
of 25.26 uM, 2.26 uM and -6.94 uM respectively (Figs. 19, 20 
and 21). The correlation coefficients were r = -0.71, r = -0.70 
and r = 0.94 respectively. From this analysis a nutrient 
concentration was computed for each calibrated satellite 
temperature (Table III). 

Sea surface temperature, nitrate and phosphate maps were 
generated by "zoom transferring" the thermal pattern in the 
IR satellite image to a navigational chart. Each isopleth 
1s the average of the nutrient or nutrient concentration on 


each side (Figs. 22, 23 and 24). 





fare ITT 


Calibration table of the satellite image for 29 October 1980 


Color pace lite Gorseeted Nitrate Phgsphate 
Graee Coumt Temperature (°C) Temperature (°C) NO, (uM) PO) (uM) 
K 100 a2 0 87.49 eo On oy 
J 101 oe? ie Zoo 0.73 
Z 102 13.34 We Ae 3236 Oe7S 
H 103 a 16.66 4.15 Ooo 
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e EOS 2203 15.44 a7 0 0.94 
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Figure 19. Nitrate versus temperature with the 
regression line for 29 October 1980. 
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Phosphate versus temperature with the 
regression line for 29 October 1980. 
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Figure 21. Nitrate versus phosphate with the 
regression line for 29 October 1980. 
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iV.  DESCUESEON 


Inferring sea surface nutrient concentrations from satellite 
derived thermal patterns appears feasible in regions of strong 
upwelling, where the nutrients to temperature correlations 
are very high. However, some errors are to be expected. The 
computer program developed by Lundell (1980) to navigate on 
the picprint can introduce an error in the position of the 
cruise track on the order of two lines and two pixels apart 
(distance ~ 2 km) [Lundell, 1980]. With the identification 
of three landmarks (Pt. Sur, Pfeiffer Pt. and Pt. Lopez) on 
Peemercou int, the author coula@” compute a ©.m.S. error of ca. 
femeines and 1.3 pixels for the 9 June 1980 cruise track, ca. 
1.0 lines and 2.1 pixels for the 11 June 1980 cruise track, 
fee Oo Pines and 1.8 pixels fox the 29 October 1980 cruise 
track. These errors are smaller than expected due to the 
proximity of the landmarks, and they will not significantly 
influence the accuracy of the thermal calibration. However, 
the process of averaging the transmittance and path radiance 
values obtained for each sharp temperature gradient along the 
PoemrmensOL the cruise track will add a significant error. 
Taking this into account, the author computed a std. dev. 
mmeeo C,e=O0.85 °C and +0.56°C for the calibrated temperature 
values on 9 June 1980, 11 June 1980 and 29 October 1980, 
respectively. From these values we can conclude that the 


temperature accuracy increases if we apply the radiative 





transfer equation in regions with sharp temperature gradients, 
as in upwelling fronts. With these errors for the thermal 
maps, it would be expected to find errors of 9% and 4% for 

the nitrate and phosphate values in the maps respectively, 

on 9 June 1980, 3% and 2% for the nitrate and phosphate values 
in the maps respectively, on 11 June 1980, and 8% and 10% for 
the nitrate and phosphate values in the maps respectively, 

on 29 October 1980. Also, there are small additional errors 
due to imperfect temperature vs. nutrients correlations, that 
are negligible when applied in regions of strong upwelling. 

The quality and the accuracy of the results depends upon the 
uniformity of the age of the upwelled surface waters (com- 

pare correlation values in June and October, in Table IV). 

It therefore seems reasonable to assume that in an upwelling 
zone, the major patterns of nutrient concentrations can be 
Zeeeeoed Using satellite IR imagery and limited in-situ thermal 
and nutrient data, and even be used in the construction of a 
biochemical model, or further, of a predictive model in frontal 


Zones. 
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VEASBLE IV 


Summary Table of Nutrients vs. 
Temperature Linear Regression 
Analysis 


DAY Linear Regression Equations Gorrce lation 


N ~ nitrate (uM) P =- phosphate (uM) 
T - Temperature (°C) 


oY June Nee —  —3eoot. + 49227 =O 
2 sae Ones /T + 4272 0). 96 
ll June Nae 4. 20 7 367 = og 7 
P= -0.29Re+ 4.35 S10) 5 S73) 


il 
I 
cme 


29 October N -27T + 2er26 -O.71 


Po =@e0.09T + Zao (06 70) 
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